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Abstract 

Recoupling matrix elements are evaluated, in the harmonic oscillator approximation, 
for all possible angular and radial excitations in processes where quarks recombine. A dia- 
grammatic representation is given. Their use is demonstrated in calculating the transition 
potential for p — > 2n in a pair creation model. 
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1 Introduction 



The description of decaying particles has already a long history. Nevertheless, only recently 
phenomenological potential models [1, 2, 3] have been proposed, for which the description of 
wave functions, hadronic widths, resonance scattering, goes beyond perturbative methods. 

These potential models are based on the 3 Po mechanism for decay. It is therefore necessary 
to have a scheme along which the coupling constants for the coupling of many quark channels 
to many hadron channels can easily be calculated. In this work we propose a scheme which is 
based on ideas which are intimately connected with the potential model of [3]. 

We assume that, during the process of the decay of a hadron, an intermediate situation 
occurs in which the original quarks and the newly formed 3 Po pair are kept together in a kind of 
hadronic soup. The formation of new hadrons out of this soup is then assumed to be caused by 
rearrangements of the quarks. 

The forces in the intermediate soup are taken to be harmonic oscillator forces. This has two 
advantages: 

• the number of possible final states is limited; 

• the calculations can easily be performed. 

Since the complications of spin are easy to handle in this scheme, because spin is not connected 
with space coordinates, we confine ourselves to the treatment of the space part of the so-called 
recoupling matrix elements 2 . 

In Sect. 2 we introduce as an example explicitly an orthogonal transformation which describes 
the rearrangement of quarks and antiquarks for the case of the decay of a meson into two mesons. 
The calculation of the recoupling matrix elements for the general case is performed in Sect. 3. The 
central result of this paper is given in formula (26) and by the related diagrammatic representation 
which is shown in Fig. 1. In Sect. 4 we show how, within this formalism, we can calculate the 
spatial dependence of the transition potential for the decay of a p meson into two pions. 



2 Spin, isospin and color degrees of freedom have in the mean time been dealt with in a separate paper: 
E. van Bevcren, Coupling constants and transition potentials for hadronic decay modes of a meson, 
Z. Phys. C 21, 291 (1984), [arXiv:hep-ph/0602246]. 
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2 The decay of a meson into two mesons 

After the formation of an extra quark- ant iquark pair out of the vacuum ( 3 Po) i n a decaying 
meson, we can treat this meson as a system of four particles. We have for instance, quarks 1 and 

3 and antiquarks 2 and 4. Suppose that 1 and 4 come from the decaying meson and 2 and 3 from 
the 3 Po pair, then we study the coupling of this system to a meson-meson final state where the 
flavor content of one meson consists of flavors 1 and 2 and of the other meson of flavors 3 and 4. 

We thus study a system of four particles with constituent masses fj,i, /12, H3 and The 
positions and the momenta of the quarks and antiquarks are given by r» and Pi (i = 1,...,4) 
respectively. On grounds of our investigations [4] it is reasonable to place the particles 1 to 4 in 
independent mass- dependent non-relativistic harmonic oscillators. The oscillator frequency u is 
universal as is noticed in [3, 4]. In the centre of mass system the motion of the four particles is 
described by the Hamiltonian: 



2 , A x 2 



' 4 \ / 4 \ 

H Pi) + E ^Wi 

u=l / \i=i / 



i=i 2fk 2 (nx + H2 + H3 + A*4) 

The Hamiltonian (1) can be expressed in terms of three independent dimensionless coordinates 
and momenta. If we want to describe the situation before decay we may choose the motion of 
4 with respect to 1, the motion of 2 with respect to 3 and the motion of the pair (3 + 2) with 
respect to the pair (1 + 4). For this we introduce dimensionless coordinates p and momenta 7r: 
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The Hamiltonian which governs the situation before decay is given by formula (1) expressed in 
terms of the coordinates and momenta (P14, ttu), (P32> 7r 32) and (p 1432 , ^1432): 



n 



lfn 2 2 2 2 2l ( o \ 

-CJ |7T M + 7T 32 + 7T U32 + p 14 + p 32 + p M32 I . (3) 



Analogously we can describe the situation after decay at best with the help of the coordinates 
and momenta (p 12 , ^12), (P34; 7r 34) and (P1234, tti234), which respectively belong to the motion of 
2 with respect to 1, of 4 with respect to 3 and of the pair (3 + 4) with respect to the pair (1 + 
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2), i.e. 



^~ f 2 2 2 2 2 2*1 

^ = 9 W T12 + 77 34 + ^ 1234 + Pl2 + P34 + Pl234/ ■ 



(4) 



It is an easy task to check that the transformation which transforms (3) into (4) is given by 



\J fa (fa + fa) (fa + fa) fa \jfa (fa + fa) (fa + ^l) ^1 



\ 



AH/U2 (/U 3 + yU 4 ) fa 



\ 



i=l 



111 (fli + ll 2 ) {fa + fa) fa \jfa {fa + fa) (fa + fa) fa 



(fa + fa) fafa^2,fa 



i=l 
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i=i 



A 
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(5) 



For the matrix a it has to be understood that 

/ a n a 12 a 13 \ ( p 14 \ 



/ P12 \ 
P34 

V P1234 / 



or 



/ 7T12 \ 1 

7T34 
V 7T1234 / , 



«21 «22 CK23 



P 3 2 



or 



/ 7T14 \ ' 

7T32 
V 7T1432 j , 



(6) 



V «31 «32 «33 / V Pl432 / 

The transformation a is orthogonal. 

For the volume element which determines the normalization of the wave functions, we choose 

,3, \ 73 / v ,3 (far?, + far 4 far 1 + fi 2 r 2 \ 

d 6 (r 2 - ri) <T (r 4 - r 3 ) <T 1 



fa + fa 



fa + fa 
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i=4 



fafafafa 



3/2 



d 3 p l2 d 3 p u 0? 3 Pl234 



Y t*i 

i=4 



fafafafa 



3/2 



d 3 p u d 3 p 32 d 3 p U32 . (7) 



The second equality in (7) follows because the transformation a is orthogonal. 

Properly normalized solutions for a given eigenvalue E of (3) now take the form 

-3/4 



^{nAm} (Pl4 ) P32 ) Pl432) 



U) 3 Y & 

i=4 
fafafafa 



H,ti,mi (Pu) 



[P32) 



^Pl432) 



(8) 

where {n,£,m} stands for the set of quantum numbers m, £1, mi, n 2 , £ 2 , m 2 , n 3 , £ 3 and m 3 , 
where ^ 

E = p " + U + , (9) 



and where the harmonic oscillator wave functions are defined by 



2r (n + £+l) 
T(n + 1)T(£ + 



1/2 



/y r ^(f)e-^ r2 1 F 1 ( 



-; r 
2' 



(10) 



Analogously, we have a normalized basis for the solutions of equation (4): 

-3/4 



X{n',l',m'} (Pl2 ; P34 ) Pl234; 



V* 

i=4 



^1^2^3^4 



V^i.S Pn^.mj, (P34) (Pn' 3 ,i' 3 ,m' 3 (Pl234; 



11^ 



The decay of a meson to a well-defined meson-meson final state, is now described by matrix 
elements which determine the decomposition of the ^-solutions (8) on the basis of X-solutions 
(11). In this way we translated part of the problem of finding a transition potential which 
describes decay, into the problem of finding decomposition matrix elements D for an orthogonal 
transformation a: 



{n,£,m} 



(Pl4 i P32 1 Pl432) 



(12) 



^ 'D{{n,e,m},{ri,e',m'}}(Pl4> P32 > Pl432 > Pl2 ' P34 ) P1234) ^{n',£',m'} (Pl2 ' P34 ) P1234) ■ 

{n',£',m'} 

Depending on the specific situation under consideration, one can always find such an orthogonal 
transformation which describes the recombination of any system of quarks and antiquaries. 
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3 The decay matrix elements 



In this section we assume that the independent motion in the centre of mass system of N + 1 
spinless particles, which move under the influence of harmonic oscillator forces, can be described 
by means of N independent and dimensionless coordinates r t and momenta pi (i = 1,...,N). 
Before the recombination, the Hamiltonian of the system reads 

H = E^P? + r?) ■ (13) 
i=i z 

Furthermore we assume that the recombination of the system can be described by an orthogonal 
transformation (repeated indices imply summation) 

r'i = aijTj and p\ = aijPj . (14) 

After the recombination the Hamiltonian (13) obtains the form 

H = E \ {ft + rl) . (15) 

8=1 

Solutions of (13) are given by 

N 

i^fnAm} (M) = IT </>mA,mi (n) , (16) 

where {n,£,m} stands for the set of quantum numbers n\, £\, mi, . . ., n^, £^, and {r} for 
r 1; . . ., r N . The relation between {n,£,m} and E reads 

N , 



E = J2 (2n i + £ i + l) . (17) 



i=l 

Analogously, solutions of (15) are of the form 



i=l 

where 
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E = £ (2n> + £> + - j . (19) 
Now, we want to generalize the X>'s of (12) to N oscillators: 

'D{{n4,m},{n'4',m'}} ({ r } i r '}) = i^{n' ,m'} ({ r '}) ^{n,l,m} ({ V }) ) ■ (20) 

Thereto we introduce the generating function of harmonic oscillator wave functions 

G(e,r) = e s 2 + 1s-r-\r* _ (21) 
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In [5] it has been shown that 

G(s,r) = £ /3(n,£) S 2 ^yr 

where 

0(n,*) = (-l) n 



(22) 



2tt 3 



1/2 



r(n + i)r(n + ^ + 

We first notice that for the orthogonal transformation (14) follows 



A? 



N 



N 



Si , r, 



i=l 



t=l 



JV Af 

t=i j=i 



(23) 



If we then multiply with the product (si) ... YJ? (sn), integrate both sides of Eq. (23) over 
the volume dQs! ■ ■ ■ dQs N and equate the coefficients of the equal powers of si, . . ., sn, then we 
find for (18) the expression 



^fn',£',m'} ({ r '}) 



(24) 
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(4tt; 
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A 
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n 



2^ + 1 



1/2 



m a • • • m i7V 



II '(n,,.!,,) ( 



m, / i=i 



Here {77.^, ^y, m^} stands for the set of 3 times N x N quantum numbers nu,4i)Wiu, . . ., 
i^nn, ^NN, m NN- An exception occurs if one (or more) of the matrix elements a^- equals zero. In 
this case it is sufficient to replace (atij) 2ni3+li:i by S^pS^^o, which replacement can be understood 
if one puts s = in (22). The symbol between brackets in (24) is defined by 



in 
ma 



(47T 



hn 
m iN 
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With [6] this can be written as a sum of products of S'0(3)-Clebsch-Gordan coefficients: 
tx ■■■ In 



mi ■ ■ ■ mN 



(25) 



in 



H+l 



N ~ 1 ( Li 



Li h+i L i+ i 




The matrix element (20) can be written in the form (using (24), (25) and the result of the 
Appendix A) 

|JV(iV-l) 



^{{n,l 1 m},{ii'/,m'}}({ r } i { r '}) 
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1=1 



3 

i i 2 



21 + 1 



1/2- 



3=1 



T{ n>+l)T[n>+£>+l 
2l> + 1 



1/2 



E 



JV IN 

11^(2^+^,2^1; 

j=i \j=i 
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(26) 



i=l j'=l J- 

This expression is the central result of this paper. 

With an extension of the rules defined in Ref. [5], formula (26) can be represented by a diagram 
as depicted in Fig. 1. Here, every left vertex {rij,£j,mj}, where the internal lines, (ny, £y, my), 
. . ., (riNj,£Nj, m Nj), end, contributes a factor 



T(n j + l)T(n j + e j + 
2£j + 1 



1/2 
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N.j 
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and every internal line, {riy, i^, contributes a factor 



7T \ 2 



O 



2£ {j + 1 



r (ray + 1) r (n 4J + t i3 + 1) 



(29) 



ni, li, mi 



\ - 

\ 

\ 

\ 



' V 

/ 

/ 

/ 

/ 



- - - - 22- -\ /--22 -^r 



/ ^--^ n 2 ,^2' m 2 



\ / 
\ \ / / 

\ v / 

N X ' 

\ * \ / 

< A 
/ \ / \ 

/ V N \ 



/ / 
/ / 
/ ' 
/ / 



\ \ 
\ \ 
\ \ 



\ 



\ 



/ / 
/ / 

.'JV. 



\ \ 
\ \ 

N S 



njy, £ N , mjy / 



^ NN NN -X 



n iV' ^AT' m AT 



Figure 1: Diagrammatic representation of formula (26). 
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4 The transition potential for the decay p — > 7T7T 

If, just before the decay of a p meson, a 3 Jo pair has been created, we are in the situation in 
which two particles (coming from the p meson) are in a relative S wave, in which the 3 Po pair 
is obviously in a relative P-wave and in which as we shall assume, both pairs are in a relative S 
wave. Just after the decay of the p meson there are two pairs (the pions). The particles inside 
each of the pions are in relative S waves, and the pions are in a relative P wave. The diagram 

(n,0,0) 
P 



(0,1, m) 
3 Po 



relative (0,0,0) 
motion 



(0,0,0) 

.__ 0,0,0——/ 7T 



/ - - 0,Q 

(0,0,0) 

0,0,0— > 7T 



'^-^s (n,l,m) relat 
0,0,0 - - - - v 



lve 
motion 



Figure 2: Diagram for p — ► irir. 



for this rearrangement is shown in Fig. 2, from which one may conclude that we only allow radial 
excitations in the p meson initial state and in the relative final-state motion of the two pions. 

Because we assume equal masses for the participating quarks, the orthogonal transformation 
(14) takes the form: 



a 



( i 

2 

1 

2 



V 



o 



(30) 



/ 



There is only one non-vanishing contribution to the sum of possibilities of the internal lines in 
the diagram of Fig. 2 (this possibility is shown!). 
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For the calculation of the diagram we need the following vertices 

(0,0,0) (0,0,0) 
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and the internal lines 



(0,0,0) 



(n,0,0) 



7T \ 2 



4/ r(n + i)r 



7? 



(0,1, m) 



V2 



With the use of formulae (31) and (32), we obtain for the diagram the result: 



l\ n 2n + 3 

2J 
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The result (33) can be used for the calculation of a transition potential. In a specific model 
we consider the quarks in the p meson as pointlike particles. The initial state is thus described 
by a wave function (cp \x), where x is the separation between the quarks. 

Similarly we treat the final-state pions as pointlike objects for which we take the same spatial 
separation as for the quark and antiquark in the initial state. The rearrangement potential which 
follows from this treatment equals 

(x\V\x') = {x\n)(n\V\n'){n'\x') . (34) 

ra,n' 

In (34) we have inserted a complete set of harmonic oscillator wave functions 

, . . -iywi 2 , / 3 m 2 \ 

{x \n) oc e 2 * <p I — n ; - ; -^wx J ■ (35) 

because the p meson is in an S wave, 

i n i\ i —\rnujx' 2 ,( i 5 /2 \ , s 

(n \x ) oc x e 2 I — n ; - ; mux J . (3d) 

because the pions are in a relative P wave. The masses in (35) and (36) are such that a pion is 
twice as heavy clS db (constituent) quark. 

If we insert in (34) for (n \ V\ n') the diagram of Fig. 2, we clearly see that n' = n. Using (33), 
(35) and (36), we obtain for (34) the result 



1\" /2n + 3 , -Imw^ + x' 2 ) ,( 3 m 2 \ / 5 
- \ x p. * \* /ml —n : — : — ujx (h —n : — 



(x\V\x') oc ^-J y — - — x'e 2 V2 ' J (j) i^- n ; _ ; — ux J ; _ ; mux > 

It is common practice to approximate this by a local one: 

V(x) = (x\V\x')5(x - x) 

in this case we obtain 



(37) 



l\ n /2n + 3 -3 muJX 2 ,( .3_m 2 \ / ^5_ 2 



V(x) oc ^-j \j—^ xe V (j){-n; -ux'j^^-n; -;mLUX 2 ) . (38) 
The form of the potential (38) is in agreement with the phenomenological potential of [3] 

V(x) oc xe^™ 2 . (39) 
where the authors quote j3 = (about) 1.7. 



12 



5 Conclusions 



The advantage of the use of harmonic oscillators for the calculation of strong decay potentials 
is not only that the calculations are very easy and that the number of possible final states is 
limited, but also that it is possible to determine the spatial dependence of the strong decay 
potential. Whether or not this spatial dependence is in agreement with experiment, might be 
checked within the decay model proposed in [1, 3]. However, since the position of the maximum 
of (38) is in accordance with the maximum of (39), not much disagreement has to be expected. 
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A Appendix 



The derivation of (26) involves the following product of integrations: 



JV 



n d 



3=1 



N 

n 

i=l 



■'n i ,e i ,m i \ r j) e 



-l(N-l) r) 



(40) 



which equals to (using Eq. 25): 

N ( / N n 

II ME (2n« + k) , 2n, + i 3 



u=l 



1 x ^(JV-i) 



T(n 3 + l)T[n J +£ J + 
2£j + 1 



Zlj ■ ■ ■ tNj 

mij ■ ■ ■ m N j 

1/2 



(41) 



m,- 



JV 



n (-i) ni 



+ 1 



2=1 



+ 1) r (n. 



Because the derivation of the result (41) is rather lengthy, we will not present its full details, but 
only outline a possible procedure. 

The delta functions in (41) follow from the fact that the functions 



e 2 



k(N-l) r 



N 

n 

i=l 



(42) 



are eigenfunctions of H 3 = | ( + p A , with eigenvalues 



N , 3 
E 2^ + ^ + - 

From the definition of the harmonic oscillator wavefunctions (10), it can be seen that the in- 
tegrations over the solid angles are exactly those already defined in Eq. (25). The remaining 
parts of the integrations may be performed by introducing the series expansions of the confluent 
hypergeometric functions. This leads to standard integrals: 



r| drj 



N 



nfo 



i=i 



3 \ _l r 2 

lj p 1 (_ n ..£. + _. r Jj e 2 , 



r + 1) 

r (n,- + i 3 - + | 



r + 1) 



i=l r {—Tlij, 



E 

fei,...,fejv 



n 



r - ra. 



i r(^ + i)r + ^ + 



JV 



JV 



r (g ^ + i%) + ^ + 5 



i=l 



i 



2r (H -E (** + ^)) 



At this stage one must also use the fact that the angular momenta are restricted because the 
€s in the symbol defined in Eq. (25), must add up to an even integer. Of the resulting product 
of sums only one term of each sum contributes, which leads finally to the expression (41). 
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